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1 Introduction 

The mean field approach when combined with effective interactions Uke 
Skyrme or Gogny is central to the understanding of nuclear structure as 
it provides the right magic numbers and a rather good description of nu- 
clear properties like masses, radii, etc, see [1] and references therein. In this 
approach nucleons are assumed to move in orbits created by a common po- 
tential and therefore the nuclear wave function for the ground state can be 
represented as a Slater determinant built upon the orbitals occupied by the 
nucleons. The common potential is determined from the effective interac- 
tion by solving the Hartree- Fock (HF) equation. On the other hand, most 
of the nuclei show, at least in their ground state, the phenomenon of nu- 
clear superconductivity due to a part of the interaction known as pairing 
interaction. When this is the case we have to introduce the concept of quasi- 
particles (given by the canonical Bogoliubov transformation) and the mean 
field wave function now becomes a product wave function of annihilation 
quasi-particle operators. The quasi-particle amplitudes are determined by 
solving the Hartree- Fock- BogoHubov (HFB) equation. 

The nuclear mean field has a strong tendency to show the phenomenon 
called "spontaneous symmetry breaking" [2, 3, 4] that appears when the HF or 
HFB wave functions do not respect the underlying symmetries of the Hamil- 
tonian. In fact this is the case, for example, for the HFB wave functions as 
they do not have a definite number of particles, i.e. they spontaneously break 
the "number of particles" symmetry. As the atomic nucleus is a finite sys- 
tem the spontaneous symmetry breaking mechanism is a mere artifact of the 
mean field approximation to generate correlations (although it allows an in- 
tuitive understanding of some nuclear structure effects like rotational bands) 
contrary to what happens in quantum field theory where it represents a real 
effect due to the infinite number of degrees of freedom. Usually in nuclear 
physics the spontaneous symmetry breaking has to do with spatial symme- 
tries like the rotational or parity symmetries. In the former case it leads to 
the concept of "deformed mean field" where the common potential felt by 
nucleons is not rotational invariant (i.e. it is deformed in opposition to a 
spherical -rotational invariant- potential) and, as a consequence, the ground 



2 J. L. Egido and L. M. Robledo 



state wave function is not an eigenstate of the angular momentum operators 
and Jz (i.e. the Casimir operators of the rotation group). However, the 
real wave function of the nucleus is an eigenstate of angular momentum and 
therefore it is necessary to go beyond the mean field approximation in order 
to have the right quantum numbers. The procedure to restore the symmetry 
is known as the Angular Momentum Projection (AMP) method [5] and relies 
in the fact that when the "deformed mean field" wave function (in the follow- 
ing the "intrinsic" wave function) is rotated the corresponding intrinsic mean 
field energy remains the same. Therefore, a suitable linear combination of 
such rotated "intrinsic states" will recover the angular momentum quantum 
numbers of the wave function and at the same time will reduce the energy 
(what is worthy from a variational point of view). Usuahy it is said in the 
literature that the "deformed mean field" wave functions belong to the "in- 
trinsic" frame of reference whereas the projected wave functions belong to 
the "laboratory" frame of reference. In order to go from the intrinsic frame 
to the laboratory one the fiuctuations in orientation have to be added to 
the intrinsic wave function in exactly the way as it is done in the Angular 
Momentum Projection framework. 

Angular Momentum Projection has been a goal of nuclear physicist for 
many years and only recently with the new computer facilities has become 
a reality for involved forces. Apart from the traditional motivations a con- 
siderable effort has been made, in the last few years, in order to implement 
Angular Momentum Projection with realistic effective interactions. The main 
reason, apart from the genuine desire of always having the best theoretical 
description, is the wealth of new experimental data in exotic regions of the 
Nuclide Chart far away from the stability line and coming from the amazing 
and very sophisticated experimental setups that have been assembled in the 
last few years. In many cases, the experimental results can not be reproduced 
or even understood adequately with a mean field description and therefore, 
in addition to considerations concerning the suitability of such interactions 
away from the stability line, effects beyond mean field have to be explored. 
As a result of the AMP calculations it turned out that away of the stability 
fine the potential energy surfaces are very soft and that shape coexistence is 
a phenomenon rather common in these exotic regions. In these cases one is 
forced to perform, besides the AMP, configuration mixing calculations. The 
most effective way to consider shape mixing is the Generator Coordinate 
Method (GCM) in which the relevant coordinates (in general the multipole 
moments) are used to generate the corresponding wave functions. The mix- 
ing coefficients are obtained by solving the Hill- Wheeler equation [6]. The 
combined AMPGCM with effective forces is a very powerful method which 
has allowed, as we will see, to understand and to predict many new features. 

Obviously there are many more methods apart from the one to be dis- 
cussed here that are successful in the description the nuclear structure phe- 
nomena at low energies. We have the traditional sheh model (see [7] for a 
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recent review and [8] for the state of the art implementations of this method) 
where first a reduced set of orbitals is chosen as to be tiie one playing the 
most important role in the physics to be described and then the full diago- 
nalization of the Hamiltonian in the space of multi particle-hole excitations 
coupled to the right quantum numbers is performed. Very good results are 
obtained with this method when the effective interaction used is fine tuned to 
the set of nuclei to be described and also when the physics to be described lies 
within the configuration space chosen. However, there are several drawbacks 
for this method: a) The dimension of the Hamiltonian matrix to be diagonal- 
ized dramatically increases with the size of the configuration space restricting 
its applicability to mass numbers smaller than 70 and also to states that do 
not involve different major shells at the same time, b) Its success depends 
upon a careful fitting of the interaction and therefore it is not suited for ex- 
ploratory calculations in new regions of the Nuclide Chart, c) Finally, it is 
difficult to recast the results in terms of traditional concepts based on a mean 
field picture of the nucleus. 

Another interesting approach is variational approach of the Tiibingen 
group [9, 10] where the lab frame wave functions of the nucleus are con- 
structed by projecting, onto all the preserved quantum numbers, an unre- 
stricted intrinsic mean field wave function. The method is sometimes com- 
plemented by allowing also multi particle-hole excitations which are also vari- 
ationally determined and then the set of wave functions obtained is used to 
diagonalize the Hamiltonian. This method fully shares one of the drawbacks 
of the shell model approach, namely the one denoted by b) in the previ- 
ous paragraph and it also partially suffers from the drawback a) of the shell 
model, namely that not so big configuration spaces can be used in the cal- 
culations. Up to now it has been possible to study nuclei in the Kr (Z=36) 
region with this method. 

An approach that can be extended to heavy nuclei is the Projected Shell 
Model (PSM) of Hara and Sun [11]. In this method the Hamiltonian is taken 
as the one of the Pairing+ Quadrupole model with single particle energies fit- 
ted to experimental data. The HFB ground state and many multiquasiparticlc 
excitations are projected onto good angular momentum and the Hamiltonian 
is diagonafized in the resulting basis. One of the advantages of the method 
over the two previous ones is that bigger configuration spaces can be used 
(up to three harmonic oscillator major shells) but it shares the deficiency 
related to the fitting procedure of the interaction in order to get a reasonable 
description of experimental data. 

Finafiy, we have the Monte Carlo Shefi Model (MCSM) approach of Ot- 
suka and collaborators [12]. In this method intrinsic wave functions are gen- 
erated stochastically and then projected onto the right quantum numbers. 
The resulting configuration is kept if its projected energy is lower than a 
given threshold generating in this way a basis of projected wave functions. 
The Hamiltonian is diagonafized in this basis at the end. As in the other 
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approaches, rather limited configuration spaces can be used (although it is 
possible to carry out calculations for nuclei as heavy as Barium ) and the 
matrix elements of the interaction have to be carefully fitted to the region of 
interest. 

The four approaches just described share the common problem of the 
effective charges that have to be introduced in the calculation of transition 
probabilities as a consequence of the limited configuration spaces. 

The main advantage of using beyond mean field approaches with effec- 
tive forces over the other approaches mentioned before relies in the universal 
character of the forces. The phenomenological effective forces are supposed 
to be valid all over the Nuclide Chart and are also supposed to contain all 
the ingredients needed to describe well low energy nuclear structure phenom- 
ena. In addition, the effectives forces are defined over the whole configuration 
space (in principle over the whole Hilbert space, including even continuum 
states) and therefore there is no need to specify which orbitals will play a 
role when a new phenomenon has to be described as all of them enter the 
game and its role will be determined by the optimization of the energy. As a 
bonus, no effective charges are needed in the calculation of transition prob- 
abilities. Another important advantage of the method comes from the fact 
that the starting point is always the mean field and therefore the results are 
much easier to interpret in terms of familiar quantities. Up to now the whole 
AMPGCM has been performed with the quadrupole moment and restricted 
to the axially symmetric case, which hints to the computational drawbacks 
of the method. In principle one would hke to take as many generator coor- 
dinates as possible, however for the nowadays computers a two dimensional 
AMPGCM sets the limits of reasonable calculations. The reason is that the 
effective forces are defined in the whole configuration space. As a consequence, 
it is capital for the method to perform reasonably well to have a good guess 
of the relevant degrees of freedom to describe a given phenomenon. Global 
properties as binding energies, quadrupole moments, etc, as well as ground 
and collective excited states are usually very well described by the method. 
However, a very accurate description of any general excited states beyond 
the AMPGCM in its present form. 

In the following we will discuss AMPGCM with effective forces of the 
Skyrme and Gogny type. In Sect. 2 the mean field approach will be briefly 
reviewed and the Skyrme and Gogny forces introduced. In Sect. 3 the techni- 
calities of AMP will be discussed and some examples will be used to illustrate 
the procedure. In Sect. 4 configuration mixing with AMP will be presented an 
its main outcomes will be discussed with an example. In Sect. 5 an account 
of the most remarkable results obtained so far with both the Skyrme and 
Gogny interactions will be prosontod and briefly discussed. In Sect. 6 we will 
end up with an outlook of the results discussed and we will discuss further 
developments of the theory. 
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2 Symmetry Breaking Mean Field 

The mean field HFB wave functions are determined in terms of quasi-particle 
creation a+ and annihilation = (a^)^ operators (Bogoliubov canonical 
transformation) 

k 

what are given as linear combinations with the amplitudes Uku and Vkij, of 
convenient creation and annihilation single particle operators c^and Ck (usu- 
ally creating or annihilating Harmonic Oscillator eigenstates). The ground 
state wave function is the product wave function defined by the condition 
a^lv?) =0 and given by 

where |0) is the true vacuum wave function and the product runs over all 
quasi-particle quantum numbers leading to a non zero wave function |(^). 
The and Vfe^ amplitudes are determined by requiring the HFB energy to 
be a minimum what leads to the well known HFB equation 

\-A* -h* ) \V U* ) \VU*)\Q-E) 

This is a non linear equation as the HFB fields h and A depend on the 
solution through the density matrix and pairing tensor (see [2, 3, 4, 1] for 
technical details on how to solve this equation). As was already mentioned 
in Sect. 1 the solution of (1) does not preserve in many cases the symme- 
tries of the Hamiltonian and is very usual to find that it breaks rotational 
invariance leading to a deformed matter density distribution which is char- 
acterized by its multipole moments like the quadrupole moments ^2/1, the 
octupole ones (73^, the hexadecapole ones 54^, etc. In order to better charac- 
terize the symmetry breaking solution it turns out to be convenient to study 
the HFB energy in the neighborhood of the self consistent minimum. In this 
way we can study whether the solution corresponds to a well developed min- 
imum or if there are other local minima around at an energy relatively close 
the one found indicating thereby that the mean field solution obtained has 
got chances of being unstable when additional correlations are included. The 
best way to do such an study is to carry out constrained calculations where 
the minimum of the HFB energy is sought but with the constraint that the 
mean values of relevant operators take a give value. In the case at hand, it is 
customary to constraint in the mean value of the mass quadrupole operator. 
The constrained HFB equation is the same as the unconstrained one except 
for the substitution h ^ h — Xo where o stands for the matrix of the single 
particle matrix elements of the constraining operator and A for the chemical 
potential that is determined as a Lagrange multiplier to force the solution to 
satisfy the imposed constraint. 
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For the interaction to be used in the HFB calculation there are several 
possible choices but nowadays the most popular options seem to be the effec- 
tive density dependent interactions of Skyrme or Gogny type. Both are very 
similar in structure the only difference being that the former is a zero range 
force contrary to the latter and therefore it cannot be used for the calculation 
of the matrix elements entering into the definition of the pairing field. 

Both the Skyrme [13, 14, 15] and Gogny [16] interactions are two body, 
density dependent phenomenological interaction given by the sum of four 
terms 

V{1, 2)=Vc+ VlS + VdD + VCoul- 

The Coulomb field t^coui = e^/ki — 'r'2\, the density dependent interaction 
VDD = ^3(1 + PcyXo)S{ri - r2)p"{{ri + r2)/2), and the two body, non- 
relativistic spin orbit potential vls = iWLsi'^i2S{ri — A Vi2)(cri + (72) 
are the same or very similar in both interactions (they usuahy differ in the 
power a entering the density dependent term). However, the central poten- 
tial vc is different in both interactions being of zero range in the case of the 
Skyrme interaction 

wc(Skyrme) = to{l + xoPa)S{ri - 

+ iti(l + xiP„)[k'l25{ri - ra) + 5{ri - r2)feL] 

whereas it is of finite range in the case of the Gogny interaction 

vc{Gogny) = ^ e-^^'-^^^" ' [W, + BiP, - HiPr - MiP^P^) . 

j=l,2 

Both contain the usual combinations of spin and isospin projector operators 
Pa and Pr- In the case of Gogny the finite range is modeled by the sum of 
two Gaussian of different ranges. The finite range was introduced to prevent 
the ultraviolet catastrophe that shows up in the evaluation of the pairing 
potential. As a consequence, it can be used for both the particle-hole and 
particle-particle channel of the HFB equations. This is not the case for the 
Skyrme interaction and extra pairing interactions with their corresponding 
parameters are usually introduced to deal with the particle-particle channel. 
Concerning the parametrizations of the force, in the case of Skyrme there are 
over a hundred parameterizations (hke Sni[17], SkM[18], SkM* [19] or the 
Lyon forces [20, 21]) that more or less yield the same bulk nuclear properties 
at the saturation density but differ in other aspects hke their behavior in 
infinite nuclear and /'or neutron matter or their performance to describe the 
collective response of the nucleus. For the Gogny interaction there are very 
few parametrizations and in most calculations only the so called DIS is used. 
It was adjusted back in 1984 [22] to reproduce nuclear matter properties as 
well as bulk properties of some selected finite nuclei. As this parameterization 
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has performed rather well in all the cases where it has been applied it has 
been kept over the years. However, taking into account the poor performance 
of Gogny for neutron matter a new parameterization was proposed [23] but it 
has not thoroughly been tested up to now. A recent review on the mean field 
properties of cfTcctivc density dependent interactions like the ones discussed 
here can be found in [1]. 

Since we are focusing in this lecture on the rotational symmetry breaking 
mechanism, which is mainly characterized by the mass quadrupole moment 
of the density distribution, the starting point of the calculations (both with 
the Skyrme or Gogny interaction) is a constrained HFB calculation with 
the mass quadrupole components as constrained quantities. As the resulting 
wave functions will be used in the context of Angular Momentum Projection, 
which represents a very tough computational problem, the calculations will 
be restricted to axially symmetric, parity conserving configurations denoted 
by \'p{q2o)) (i.e. by construction only the /x = component of the quadrupole 
tensor can be different from zero). These HFB wave functions are obtained 
as a solution of the corresponding HFB equation with the constraint in the 
mean value of the /it = mass quadrupole operator (</5(g2o)| -2^ — 1/2(2^^ + 
y^) \^{Q2o)) = 920- Also in both kind of calculations it is customary to neglect 
the contributions of tlic Coulomb field to the exchange and pairing potentials 
due to the high computational cost associated with the calculation of those 
fields (see [24] for a thorough discussion of these contributions) . In the case 
of tlie Gogny force calculations we have added the exchange Coulomb energy 
evaluated in the Slater approximation at the end of the calculation in a 
perturbative fashion. 

Usually the HFB equation is solved for the Skyrme interaction in coor- 
dinate space representation by introducing a Cartesian mesh whereas in the 
case of the Gogny force a Harmonic Oscillator (HO) basis is used including 
up to A^o complete shells. The HO length parameters are chosen to be equal 
in order to preserve the rotational invariance of the basis which is a very 
important requirement for the subsequent Angular Momentum Projection 
calculations [25]. 

As an example of the kind of results it is possible to obtain we have plotted 
in Fig. 1 the HFB energies computed with the Gogny force as a function of 
q = q2a/A^'''^ (the variable q is defined to be independent of mass number 
A, like the deformation parameter /3) for five representative nuclei, namely 
^^Mg, ^^S, ^^^Er, ^^^Pb and ^°^Pb. These nuclei have been chosen as they 
represent typical cases: In the doubly magic ^°*Pb there is no spontaneous 
symmetry breaking and the energy has a minimum aX q = with a very deep 
and stiff weh. In the nucleus ^^^Er we have two minima one prolate {q > 0) 
which is the absolute minimum and the other oblate {q < 0) separated by a 
rather high barrier. In the nucleus "^^S we have an spherical minimum but at 
a given deformation an excited local minimum appears corresponding to a 
super-deformed configuration (obviously there are better examples of super- 
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Fig. 1. HFB energies as a function of the quadrupole deformation parameter q = 
920/^^''^ for five relevant nuclei. Energies are relative to the respective absolute 
minimum 

deformed statos both in tho rare earth region as well as in the Hg region; 
we have chosen this particular example as is the only one where subsequent 
AMP calculations have been carried out) . In the case of the neutron deficient 
magic "'^^^Pb we have three minima one spherical (the absolute minimum) and 
two excited minima, one prolate and the other oblate, which have energies 
very close to the one of the ground state (shape coexistence) . Finally, in the 
neutron rich iV = 20 magic nucleus ^^Mg we obtain a spherical minimum 
{q = 0) and a shoulder at around 1 MeV excitation energy. Similar results 
are obtained for the Skyrme interaction as we wih see in Sect. 5. 

3 Angular Momentum Projection 

In the following subsections we wih introduce all the technology to carry out 
angular momentum projection without entering into much details as most 
of them can be found in the literature. However, we will devote some more 
attention to the issue of how to deal with density dependent interactions. 
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3.1 The Projector Operator 

Wave functions | ^im) eigenstates of and Jz with eigenvalues + 1) 

and KM respectively can be built out of a given deformed mean field state 
I <p) by applying the angular momentum projector 

\^im)=Y,9'kPLk\v)- (2) 

K 

The angular momentum projection operator -P^^ given by [5] 

PLk = j dnv'^MR{Q) (3) 

where f2 represents the set of the three Euler angles {a,(3,j), 'D\^j^{Q) are 
the well known Wigner functions [26] and R{fi) = er^'^''' er'^^'^y er^^'"'' is the 
rotation operator. That | "Z'/m) is an angular momentum eigenstate is very 
easy to check. Applying the rotation operator R{fi) to it we obtain 

R{Q) I ^im) = '^^Y.9k j dn'V'j^j,{n')R{n + n') I ^) . 

Now the integral on the three Euler angles is rewritten as 
J dQ'V'^Ki^' - f2)R{n') I ^) = j dn'Vi\{Q')R{Q') \ ^) 

which shows that 

R{^2)\^IM) = Y,V'MLm\'^IL) 
L 

The quantities g]^ are arbitrary at this point and are usually determined in 
a variational sense to yield a minimum of the projected energy. The resulting 
equation will be discussed in the next subsection. Other relevant properties 
of the projection operator can be easily deduced from the following represen- 
tation of the projector [27] 

a 

Using this representation is very easy to derive the property 

[Pmk)^ Pm'K' = ^II'^MM'PkK' (5) 



as well as 
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{Pmk) — Pkm ■ (6) 

An interesting quantity is the overlap between the intrinsic wave function 
\ip) and the projected wave functions | •f'/Af) that gives the probabihty am- 
phtude of finding the angular momentum / and third component M in the 
intrinsic wave function. In Fig. 2 we have plotted such quantities as a func- 
tion of the quadrupole deformation for the axially symmetric HFB solutions 
obtained for the nuclei ^^Mg and ^^"^Er that were discussed in Sect. 2. As the 
intrinsic wave functions are axially symmetric the only meaningful quantity 
in this case is A/'^(g2o) = {'P{l2o)\PmW{l2o)) that is the quantity plotted. 
We observe that the spherical solution (520 = 0) only has / = component 
or in other words the HFB solution corresponding to 520 = is already an 
eigenstate of angular momentum with / = 0. For increasing deformation the 
/ = probability decreases whereas the other components increase up to a 
given point and from there on they decrease. Finally, for large deformations 
all the components are of the same order of magnitude implying a strong 
spreading of angular momentum for the intrinsic state. 




Fig. 2. Projected norms Af'{q} for 7 = 0, 2, 4, 6 and 8 (dashed lines correspond 
to 7 = 2 and 6 ) plotted as a function of the quadrupole deformation parameter 
q = q2o/A^/^ for *^Mg and ^^^Er 

Finally, it has to be mentioned that there are other representations of the 
projector P^u^ as the one of [28] that use the rotation operator in the form 
R{k) = exp{—iKnJ). 

3.2 The Projected Energy 

The projected energy is simply given by 
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and the fact that this quantity is independent of M is a direct consequence 
of the invariance of the Hamiltonian under rotations. 

Up to now nothing has been said about the gj^ coefficients as the geomet- 
rical properties of | ^im) do not depend on them. However, we can exploit 
this degree of freedom in a variational sense to minimize the projected energy 
with respect to those parameters. The resulting equation reads 

J29k' {{<fi\HPkM-E'{,p\P^^,\,p)) =0. 

K' 

To compute {(p\HP^j^, \ip) we have to perform a three dimensional integral 
over the Euler angles of an integrand that involves the well known Wigner 
functions and the Hamiltonian overlap {(f\H R{f2)\ip) . The latter quantity can 
be straightforwardly evaluated by means of the extended Wick theorem [29] 
and the sign problem of the norm dealt with as suggested in [30], but it turns 
out that its evaluation for the case at hand is a very intensive computational 
task that has to be repeated for all the Euler angle mesh points used to 
evaluate the three dimensional integral numerically. For present single CPU 
computers this is a task that demands of the order of tens of hours to com- 
plete. However, for the Skyrme interaction there has been an attempt [31] 
to carry out the full calculation for low spins by restricting the number of 
mesh points in the a and 7 integrations and this is a path worth to be ex- 
plored in the future. For the other approaches mentioned in the introduction 
the computational cost is much smaller due to the restricted configuration 
spaces and therefore the full projection is routinely carried out. Fortunately, 
in many cases the restriction to axialljf sj^mmetric intrinsic wave functions 
seems to be a sound approximation [32] that reduces the computational bur- 
den by almost two orders of magnitude. For even-even nuclei when \ip) is 
axially symmetric it satisfies Jz\(p) = and therefore {f\HR{n)\Lp) reduces 
to {ip\He~'^^'^^ \(p) which is independent of a and 7. As a consequence the in- 
tegrals in a and 7 are trivial and yields (27r)^(5if,o<^if',o- With this restriction 
we get for the projected energy 

_ d(3sm{miomv\He-^^'^\v) _ J,^ dp sm{P)dUP)h{P)n{P) 
J^dPsin{0)dU0){v\e-'f'-'yW) /;d/3sin(/3)rfo^o(/3)«(/3) 

with n(/3) = {f\e^'^^'^^\ip) and h{f3) = {(p\He^^^'-'^\Lp). Using as self-consistent 
symmetry for the intrinsic wave function the signature (essentially the inver- 
sion of the X axis, and given by the operator iTe^'"^'^"') it is possible to reduce 
the integration interval to [0, 7r/2] and to show that for even-even nuclei the 
projected energy is only defined for even values of the angular momentum 
/ when the intrinsic state is reflection symmetric. When the intrinsic state 
is not reflection symmetric the reduction of the interval of integration im- 
plies a projection onto good parity tt in addition to the angular momentum 
projection and the rule (—1)^ = tt is obtained. 
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3.3 How to Deal with the Density Dependent Term of the 
Interactions 

When dealing with density dependent (DD) forces we face the problem 
that the DD term is only defined at the mean field level: In the evalua- 
tion of the expectation value of the Hamiltonian {ip\H\ip) the DD inter- 
action is given in terms of p{R) = {(p\p(R)\'^) which is evaluated as an 
expectation value with the same wave function used in the evaluation of 
the energy. However, in the calculation of the projected energies we need 
the Hamiltonian overlap {(p\H R{f2)\(p) and, as the DD term has a phe- 
nomenological origin, it is not clear which density should be used in the 
DD term of the interaction for the evaluation of the above mentioned over- 
lap. The same problem arises with the overlap needed in the Generator Co- 
ordinate Method (GCM), namely {(po\H\(pi). Several prescriptions have 
been proposed to cope with this problem: in one [33], inspired by a general- 
ized Brueckner expansion, the density is replaced by the linear combination 
^{{ip\p{R)\ip) + {ip\f){R)R{n)\Lp)); in other called the "mixed density" pre- 
scription the quantity {(p\p(R)R{f2)\ip) / {ip\R{f2)\(p) is proposed [34] (see also 
[35] for a thorough discussion in the context of Particle Number Projection) 
and it is inspired by the way the Hamiltonian overlap is computed with the 
aid of the extended Wick theorem. There are some more prescriptions, like 
the one proposing a rotational invariant density dependent term, but we will 
not comment them here. In any case, as in most of the prescriptions the DD 
term breaks rotational invariance and even hermiticity, it is very important 
to show that the prescriptions are consistent with general properties of the 
projected energy like being a real quantity and that the projected energy is 
independent of M, i.e., Eq. (7) is satisfied. In [36] has been shown that the 
"mixed density" prescription {(p\p{R)R{n)\<fi) / {<p\R{n)\ip) satisfies the two 
previous requirements as it also does the other prescription. At this point it 
might seem paradoxical to have a Hamiltonian which is not rotational invari- 
ant but this apparent paradox can be solved if density dependent Hamiltonian 
are thought not as genuine Hamiltonian but rather as devices (in the spirit 
of the Density Functional Theory) to get an elaborated energy functional of 
the density. 

3.4 Transition Probabilities and Spectroscopic Factors 

Transition probabilities are the physical quantities that have most sensitivity 
to the approximations made to the wave functions. For instance, transition 
probabilities have selection rules that can not be reproduced unless the wave 
functions used arc cigenstatcs of the angular momentum Casimir operators 
and Jz- In this section we will present the way they are computed when 
the wave functions used come from an angular momentum projected intrinsic 
state and speciahzed formulas for the case of an axially symmetric intrinsic 
state will be presented. 
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The only property we will need is the transformation law of the multipole 
operators Qaju under rotations 

m)Q^^R\Q) = 2^mV(^)^a/ (9) 

Using the well known result for the product of two Wigner functions [26] as 
well as the definition of (3) for the angular momentum projection operator 
and the property 

PmkPr'M' = ^ij'^k,k'Pmm' (10) 
we obtain after some algebra the result 

X ^(-)'^'-m(j.^Am I IfKf)Q^^,Pl'^^ . (11) 

With the definition of the projected wave functions (2) and the previous 
result we obtain 

/* in I* \ (JiMihl^UiMll/T iin II n 

{y^ifMf I QxiM I ^i.Mi) = ^^j^ _^ ^ — (-'/ WQxlUt) (12) 

with 



^ m' -Kf) I d^'^i'li^X^f I Qx^'Ri^) I <^^)-(i3) 

The previous derivation only uses the tensor properties of the electric multi- 
pole operator and therefore also applies to the case of the magnetic multipole 
operators. 

Taking advantage of the axial symmetry of the intrinsic wave function 
as well as the self-consistent symmetry ne~™'^'' we can simplify the above 
expressions as follows. First we have 

{^f I Qv^(^) I = ^'"'"'('Pf I Qx^.'e-'^''' I ^i) (14) 

that leads to 

X /; d/3sm{P)d':\^{P){^f I Q^^'e-^f'jy \ ^f) . (15) 
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Applying this result to the expression (13) we obtain 

= (2/. + l)(2Z, + l)(-,'.-i±t)i5:(/.,A,/,) 

X rf/3sin(/3)/';,p(/3)(^/ I Q^^>e-^^-'^ \ (16) 

where we have reduced, in the last line, the integration interval to half the 
original one making use of the self-consistent signature symmetry of the in- 
trinsic wave function. 

3.5 Variation Before and After Projection 

Now that we know how to compute the projected quantities we can take for 
a given nucleus its mean field ground state and project it out to obtain the 
projected energies for different angular momenta. This procedure is called 
the Projection After Variation (PAV) method as the intrinsic state is deter- 
mined at the mean field level and projected afterwards. For the nuclei used 
as examples in Sect. 1 we will obtain for the / = ground state the same 
projected energy as the intrinsic one for all the examples discussed except 
for ^^^Er where a lowering of 3.07 MeV is obtained. For 1 = 2 the projected 
energy for those nuclei with an spherical intrinsic state is an indeterminacy 
of the type zero divided by zero. For the deformed ^*'^Er we will obtain for 
/ = 2 an excited state which is 139.4 keV higher than the / = projected 
ground state. 

As the energy gain due to projection for the ground state of ^^^Er was 
of the order of a few MeV one may wonder whether this energy gain can be 
strong enough for some deformed configurations as to overcome the mean 
field energy differences and produce projected energies which could be lower 
than the one given by the PAV method. To elucidate this possibility we have 
plotted in Fig. 3 the 7 = and 7 = 2 projected energies along with the HFB 
results (dashed line) for all the nuclei considered in Sect. 2. 

For the well deformed ^^^Er nucleus the effect of projection is just an 
overall shift of the energy with respect to the HFB curve except around 
(7 = where for 7 = the projected energy has to coincide with the intrinsic 
one. In this nucleus the 7 = 2 curve is so close to the 7 = one that is 
almost indistinguishable from the latter with the scale used for the plot, 
except around q = where the intrinsic wave function has an overlap zero 
with the 7 = 2 projected state (see Fig. 2 ). 
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Fig. 3. Angular momentum projected energies for / = and 1 = 2 along with the 
HFB energy (dashed curves) plotted as a function of the quadrupole deformation 
parameter q = q2o/A^^^ for five relevant nuclei. 

Not surprisingly the projected energies show minima in all nuclei consid- 
ered (except ^^''Er) which correspond to deformed intrinsic configurations ! 
In addition, the intrinsic configurations corresponding to the / = and 7 = 2 
minima are in the examples presented very close but they do not share exactly 
the same intrinsic state. If we now take as the wave function of the system 
for each angular momentum / the one that gives the lower projected energy 
as a function of q we will be doing things properly (the variational princi- 
ple again) and we will be using a restricted version of the Projection Before 
Variation method (PBV). In the PBV method the intrinsic wave function 
for each angular momentum is determined by minimizing not the intrinsic 
energy but the projected one. The example mentioned above is a restricted 
version of the PBV method as we are not exploring, by considering only wave 
functions constrained to a given quadrupole deformation, the whole Hilbert 
space. However, we can argue that as the quadrupole moment is the main 
quantity characterizing the amount of rotational symmetry breaking, it is 
the most relevant degree of freedom in terms of the PBV method and there- 
fore just considering it gives results very close to the ones of the exact PBV 
procedure. 
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At this point it has to be said that except the shell model, all the tech- 
niques mentioned in the introduction aiming to a better description of the nu- 
clear many body problem are different versions (or views) of the PBV method. 
What is different among them is the way the Hilbert space is searched (fully 
in the Tubingen approach, by means of multiquasiparticle excited states in 
the PSM, stochastically in the MCSM, etc). 

Coming back to the discussion of Fig. 3 it is amazing the large effect 
caused by projection for the 1 — ground state even in spherical nuclei ! 
The two light nuclei "^^Mg and "^^S become well deformed when projection is 
considered (in good agreement with the experimental results for ^^Mg). Even 
the double magic nucleus 2"Spb becomes slightly deformed^ and in this way it 
gains 1.8 MeV in correlation energy. This is an important fact (also observed 
in other doubly magic nuclei like ^^O, *°Ca and ^^Ca) that should be taken 
into account in the fitting procedure of the effective interactions treated in 
this paper. In this respect, the effect can also be relevant for the evaluation of 
masses with astrophysical purposes: in a recent fit of the Skyrme interaction 
[37] to the known nuclear masses the rotational energy correction is included 
(in an approximate way to be discussed below) for deformed nuclei in the 
spirit of the PAV and therefore is disregarded for spherical nuclei. However, 
taking into account the results just discussed for ^"^Pb the rotational energy 
correction should also be considered for spherical nuclei although the proce- 
dure to evaluation this quantity can be much more cumbersome than for the 
situation of a well deformed ground state. 

The important quantity in the description of the ground state in the con- 
text of the PBV method is the ground state energy gain due to the projection 
and given for even-even nuclei by 

Erec = Ehfb — E^~'^ (17) 

It is also called the Rotational Energy Correction (REC) and is a quantity 
which increases with increasing deformation of the system and its typical 
values are in the range of a few MeV. This is an important correction when 
shape coexistence is present in the system as it can substantially modify 
the picture emerging from the mean field and make an excited minimum or 
shoulder the ground state. It also reduces in a few MeV the fission barrier 
heights with respect to the mean field result and can have important effects 
even for spherical nuclei. Several examples will be discussed below in Fig. 6. 

^ In fact there are two deformed minima, one prolate and the other oblate with 
approximately the same energy and the same absolute value of the quadrupole 
moment. In these cases as we will see in Section 4 one should perform config- 
uration mixing calculations. The resulting wave function provides an spherical 
nucleus and an additional energy lowering of about 1 MeV 
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3.6 An Approximate Evaluation of the Projected Energies: The 
Strong Deformation Limit, the Kamlah Expansion and the Like 

The exact expressions for the projected energy and transition probabilities 
are too complicated to have an intuitive understanding of their effects. There- 
fore it is convenient to have estimations given in terms of simple quantities 
expressible in terms of the underlying intrinsic state. On the other hand, 
the evaluation of the exact projected quantities is rather costly in terms of 
computational time an is desirable to have good estimations that are eas- 
ier to compute. Such estimations can be easily obtained by means of the 
Kamlah expansion [38, 39] which is valid in the strong deformation limit. To 
understand the essence of the approximation let us consider for simplicity the 
axially symmetric case and angular momentum zero. The projected energy 
is given in this situation by 

^,^0 ^ /;/^rf/3sin(/?)fe(/3)n(/3) 
/;/'d/3sin(/3)n(/3) 

where, as compared to (8) we have reduced the integration interval to 7r/2 
by using the signature as a symmetry of the intrinsic wave function. If the 
intrinsic wave function is strongly deformed, the overlap n(/3) will only be 
significantly different from zero when /? » and will experience a fast decay 
away from that value. In fact, n(/3) can be approximated by [3] 

n(/3)«e-5<^')/''. (18) 

Next, we assume that h{(3) is a smooth function of (3 that can be expanded 
to a good accuracy in a power series up to second order in (3 

h{(3) ^ho~ ^/i2/3' + . . . (19) 
with ho = {H) and /12 = — '^ji^i^ |/3-o' '^'^^^ projected energy is given by 

I A n j2\\ fn''''^d/3sin(/3)/3=^e-5<'^y>'^^ A iU- -C ^ l 

where Aq[LJ,,)) = — -r- — -^-^ — 1,,2,,,2 • As this mnction goes to one when 

(Jy) >> 1 we finally arrive to the well known formula for the rotational 
energy correction in the strong deformation limit 

E'=° = ho-^h2 (20) 

which is nowadays widely used in many calculations to get an estimation of 
the effect of angular momentum projection in the / = ground state. To test 
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the accuracy of this approximation we have computed for the nucleus ° Er 
the exact quantities n(/3) and for different quadrupole deformations and 
the results are presented in Fig. 4 along with the approximate estimates of 
(18) and (19). 




Fig. 4. A comparison of the exact n(/3) and h{l3) (full line with symbols) with the 
approximated expressions of (21) and (22) (full line) as well as the standard ones 
of (18)and (19) (dashed lines) for different deformation parameters q. The quantity 
h{/3) has been multiplied in each panel by the indicated factor 



We observe how for strong and moderate deformations the approximated 
estimates do rather well and we can conclude that the approximation to the 
/ = projected energy is reasonable. However, for small deformations this is 
not the case and we expect a failure of the method. As it was suggested by 
Reinhard [40, 41, 32] the behavior of the overlaps can be better approximated 
for the small deformation case by using the ansatz 

n(/3) f«e-5<-^S>«''^'W (21) 

and 

h{P)-ho-^h2sm\(}) + ... (22) 

instead of (18) and (19). This ansatz is inspired by the properties n(/3) = 
n(7r — /?) and /?,(/?) = hin — /?) that the exact quantities have to satisfy when 
the simplex symmetry is imposed in the intrinsic wave function and also by 
analytical results obtained with harmonic oscillator wave functions. To check 
the quality of the new ansatz we have plotted in Fig. 4 the new quantities. 
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5 10 15 20 25 30 -0.2 -0.1 0.1 0.2 

X q{fm^) 

Fig. 5. Left panel: The function A{x) defined in (24) as a function of x. Right 
panel: The expectation value {Jy) as a function of the deformation parameter q 



We observe how the quality of the ansatz is excellent for small deformations 
and it also improves the results at modoratc and strong deformations. 
Using the new ansatz, the projected energy for / = is now given by 

E'=' = ho- lh2Ai{J^)) (23) 



with 

/;/'d/3sin(/3)/?^e- ^<-^g^ si„^(/3) 
/;/'d/3sin(/3)e- 

To understand the cliangos induced by the new ansatz the universal function 
A{x) has been plotted in the left panel of Fig. 5 as a function of x. There we 
observe how for big values of x it goes to 1 recovering the standard results 
and it rapidly decreases for x going to zero in good agreement with the fact 
that (Jy) = corresponds (for axially symmetric systems) to an spherical 
wave function where the effect of angular momentum projection is null. In- 
terestingly, the function A{x) has a maximum at a; « 7 thus enhancing the 
rotational energy correction at this point. From the shape of this curve it is 
clear that A{x) induces strong changes in the energy for x values smaller than 
7 and an smooth modulation of it at larger values. The interesting question is 
to know the deformation range for which x, i.e., (Jy), satisfies < (Jy) < 7. 
In the right panel of Fig. 5 we have represented {Jy) versus q for the indicated 
nuclei. We find that the relevant intervals are —0.2 < q < 0.2 for ^^Mg and 
-0.05 <q< 0.05 for i64Er and ^ospb. 

In Fig. 6 we now compare the exact rotational energy correction (17) at 
zero spin with the standard approximation (20) and the one by Rcinliard (23) 
for the mentioned nuclei. The exact RFC is zero for spherical intrinsic states, 
it typically increases rather abruptly for small deformations and at some point 
it slides down and from there on its growing rate stabilizes to a smaller value. 
It is interesting to see that the largest changes in the REC are given in the 
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predicted range -0.2 < g < 0.2 for ^^Mg and -0.05 <q< 0.05 for ^''^Er and 
^°*Pb. As can be observed in this plot the standard approximate expression 
of (20) does pretty bad in the corresponding deformation ranges indicated 
above while Reinhard's one does pretty well even in the critical regions. This 
is an encouraging result as it can enormously simplify the estimation of the 
rotational energy correction in all kind of nuclei, spherical or deformed, and 
therefore can make possible its use in massive mass evaluations Hke the ones 
needed for astrophysical purposes. 




Fig. 6. Rotational energy corrections as a function of the quadrupole deformation 
parameter q = q2o/A^^^ for ^^Mg, ^'^''Er and ^"^^Pb. The full line represents the 
exact result of (17) whereas the dashed line stands for the approximation (23) and 
the dotted one for the standard correction, see (20) 



What we have said up to now is based in the assumption that the intrinsic 
wave function is axially symmetric and we have restricted the discussion to 
zero spin. In the general case things can be worked out for the full projection 
[38, 39] and the following expression is obtained in the strong deformation 
limit 

+ (25) 

where Jy is the Yoccoz moment of inertia given by 
1 {AHAJ^) 



2Jy 2((ZiJJ)2 + (j2)2 + (j2)2) 

and 



(26) 



. = i^. (27) 

In order to pin down the physics behind this expression let us assume that the 
intrinsic wave function is time reversal invariant. In this case u — (Jx) = 0. 
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If, additionally, we impose an axially symmetric wave function then ( j|} = 
and the approximate projected energy becomes 

which corresponds to the energy of a rotor with intrinsic energy {H) — \jj , 
angular momentum / and moment of inertia Jy ■ In other words, we have 
a rotational band with the typical /(/ + 1) behavior and with the Yoccoz 
moment of inertia. The intrinsic energy is not simply (if) but it is reduced 

by the so called rotational energy correction (REC) given by 2jJ ■ Unfortu- 
nately, in this case the ansatz proposed by Reinhard is not easy to generahze 
and work is needed in order to get equivalent results in the most general case. 
In this respected, perhaps the representation of [28] can prove to be useful. 

Another interesting aspect of the approximate energy obtained in (25) 
concerns its use in the context of the Projection Before Variation method. 
If one assumes that the approximate energy is a sound approximation to 
the projected energy then, requiring the intrinsic state to lead to a mini- 
mum of that expression, is an approximate PBV minimization process. The 
minimization process gets rather cumbersome due to the rotational energy 
correction term but it has been shown in many examples that this correction 
stays rather constant in the domain of validity of (25) and therefore it is 
safe to neglect it in the variational process. Neglecting the variation of the 
REC one ends up with the variational equation of the self-consistent cranking 
model 

b{^p\B -u:J:c\^p) = ^ (28) 
where w is determined as to satisfy the constraint 



(<p|j.l^) = V/(/ + i)-(J|>- 

This equation has to be solved for each value of I and the intrinsic wave 
function is therefore / dependent. Another interesting aspect is related to 
the fact that Jx is odd under time reversal and therefore the intrinsic wave 
function breaks time reversal invariance as well as axial symmetry. The energy 
spectrum obtained is well represented by the energy of a rotor 

2J'scc(/) 

where the self-consistent cranking moment of inertia eventually depends on 
/. If we now compare this spectrum with the one given by (25) corresponding 
to the PAV method we notice that the two moments of inertia are formally 
different. In the cranking model we have the self-consistent moment of inertia 
whereas in the PAV method we have the Yoccoz one. How different they are 
depends on the interaction (probably on the effective mass but there is no 
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formal proof) but to give a taste of the kind of differences one can expect 
we can say that in the case of -^^"^Er the energy difference E^^"^ — E^^^ is 
139.4 keV for the exact PAV whereas the cranking model gives 78.8 keV, 
that is the cranking result is roughly a factor 0.6 smaller than the PAV one. 
Formally, the cranking results are better than the PAV (even though they can 
compare worse with experiment due to the interaction used) as the former 
are based on the PBV method. However, this does not mean that we have 
also to use the self-consistent cranking moment of inertia (or its close cousin, 
the Thouless Valatin moment of inertia) in the evaluation of the rotational 

energy correction Erec = 2jJ '^^ it is sometimes done: the REC has always 
to be evaluated with the Yoccoz moment of inertia. The above discussion is 
also in good agreement with the results of [42, 43] where also an approximate 
PBV was carried out but in a different context. 



4 Configuration mixing 

In the previous sections we have had the opportunity to glance to several An- 
gular Momentum Projected Energy Surfaces (AMPES) belonging to different 
nuclei. In many of them we have obtained several minima which are close in 
energy and are also separated by low and narrow barriers and therefore it 
is expected that configurations mixing of those states will lead to a further 
reduction of the energy of the states. In this section we will investigate this 
possibility in the framework of the Angular Momentum Projected Genera- 
tor Coordinate Method (AMPGCM). A general AMPGCM wave function 
I '^im((^)) is written as a linear combination of the "generating" functions 
I !^/M(q)) which are obtained by projecting onto good angular momentum 
a set of intrinsic wave functions which are characterized by several param- 
eter (usually multipole moments of the mass distribution ) q = {qi,qi., ...}. 
Explicitly we have 

I ^im{<j)) = [ dq/'"^(q) I ^iM{<i)) =J2[ '^q7^"(q)^MKlv(q)) • (29) 
J p. J 

The amplitudes /^''^(q) are solutions of the Hill-Wheller (HW) equation [6] 

E / rfq'7x'(q')((<^(q)|if^KK'lv(q'))-i^'"^(<p(q)l^xK'l¥'(q'))) =0. 

K' ^ 

This equation is derived by imposing the AMPGCM energy to be a minimum. 
As it was mentioned before the calculations are usually restricted to axially 
symmetric {K = 0) configurations and in this case the angular momentum 
projected Hill- Wheeler equation reduces to 

I (iq7''^(q')((^(q)l^^ool'/'(q'))-i^''"(^(q)l^ool^(q')))=o. 
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This is an integral equation where the Hamiltonian kernel H^(q, q') = 
(i^(q)|i?Pjo|(y5(q')) and the norm overlap A/'^(q, q') = (93(q)|Pjol'^(q')) play 
the central role. The HW equation is a non-orthogonal eigenvalue equation 
that is usually recasted in terms of orthogonal quantities by diagonalizing 
first the norm overlap 

j dq'A/-^(q,q')zii(q')="f^*i(q) 
what allows to write the norm overlap as 

I 

Inserting this expression into the HW equation and defining the amplitudes 
and using the inverse relation 

/'"(q') = E(;^<(q') (30) 

we end up with the HW equation written as 

Y,{Hg-E'''^5w)glr = Q 
I' 

where the collective image of the Hamiltonian kernel Hi(, is defined as 

- J ^q^q kF^^ ^'^''^ ^ ■ 

The solution of this reduced equation produces the eigenvalues E^'"' labeled 
by the cr index and the eigenvectors g^'"^ from which we can compute the 
amplitudes f^'"{c{) using (30). 

In order to interpret the GCM results it is customary to introduce the 
collective wave functions 

/■'^(q) = / rfq' (aTV^V (q,q')/^"^(q') = ^ff/x^nKq) 
■' I 

and the reason is that they are orthogonal, contrary to their counterpart 
/^''^(q')) and therefore their modulus squared can be interpreted as a prob- 
ability amplitude. 

In our particular case the set of labels q reduces to the quadrupole moment 
920- If we take into account the axial symmetry imposed to the intrinsic wave 
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functions and write explicitly the projector operator we end up with the 
following expressions for the norm overlap 

A/''(g2o,4) = (2/+ 1) /' d/3sin(/?)4*(/3)(y,((?2o) | ^'^'^ \ ^{q'20)) (31) 
Jo 

and the Hamiltonian kernel 



W'(?20,4) = (2/+ 1) /' rf/3sin(/?)4S(/3)(^feo) | HoDe-'^'^ \ ^{q^o)) 
Jo 

(32) 

where Hdd now depends on the density p^^'^{r) given by 
-GCM/^^ ^ ((^(920) I p{R)e-^P^y I ^(g^o)) 

This is the obvious generalization in the framework of the configuration mix- 
ing calculation of the mixed density prescription described in Sect. 3.3 for 
the density dependent part of the interaction. 

Finally, from the knowledge of the amplitudes /^''^(g2o), we can compute 
the reduced B{E2) transition probabilities and the spectroscopic quadrupole 
moments Q^'^'^^ {I , <j) . This is one of the main motivations for carrying out a 
configuration mixing calculation of angular momentum projected wave func- 
tions in the case of Gogny and Skyrme forces, since both interactions allow 
the use of full configuration spaces and then one is able to compute tran- 
sition probabilities and spectroscopic quadrupole moments without effective 
charges. In the framework of the AMPGCM the B{E2) transition probability 
between the states {Ii,ai) and {If,af) is expressed as (see Sect. 3.4 for more 
details) 

B{E2,hai^Ifaf) = ^^ (34) 
X J dq,dqfr'f''^f{qf){Ifqf \\ Q2 \\ hqi) f""' {qi) 
and the spectroscopic quadrupole moment for the state (/ > 2, a) is given by 



Q--(J,a)=e^/ifl(j2/j (35) 



j dqidqfr'^%qf){Iqf \\ Q2 \\ Iq^f'^qi). 
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5 Results 

In order to illustrate the whole procedure we will discuss now the ^^Mg 
nucleus in detail and at the end we will give an overview of the available 
results obtained with the Skyrme and Gogny interactions. The nucleus '^^Mg 
is a nice example as it is one of the neutron rich magic (N=20) nuclei where 
the shell closure is broken. The most convincing experimental evidence for 
a deformed ground state in the N=20 isotopes is found in the ^^Mg nucleus 
where both the excitation energies of the lowest lying 2+ [44] and 4+ [45, 
46] states and the B{E2,0~^ —>■ 2+) transition probability [47] have been 
measured. The low excitation energy of the 2+ state, the high value of the 
B{E2) transition probability and also the ratio E{4:f )/E{2f) = 2.6 are fairly 
compatible with the expectations for a rotational band. 




-1 1 2 -1 1 2 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 

Cf20 (b) Q20 (b) ^20 



Fig. 7. On the left hand side, the square of the "collective wave functions" (in 
fact, E^''^ + 25\r/'"'{q2tt)\^) for cr = 1 and 2 (thick lines) are plotted along with the 
corresponding projected energy (thin line) as a function of ^20. On the right hand 
side, the projected energies E^'"^ along with the / = projected energy are plotted 
as a function of 520. The projected energies have been placed according to their 
average intrinsic quadrupole deformation (920)^ 

We have carried out Angular Momentum Projected configuration mix- 
ing calculations along the lines described in Sects. 2, 3 and 4. The AMP 
results have been already discussed in subsect. 3.5 and the AMPGCM ones 
are summarized in Fig. 7. On the left hand side of the figure we have plot- 
ted, as a function of (720 and for / = 0, 2, 4 and 6, and for practical reasons 
a quantity related to the square of the "collective wave functions", namely 
E^ + 25|5^''^ ((720)1^, for a = 1 and 2 as well as the projected energy curve. 
By looking at the tails of these quantities we can compare the AMPGCM 
energies E^ with the corresponding projected energy and have a feeling of the 
correlation energy gained by configuration mixing. For / = we observe how 
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the ground state has an important mixing between the prolate and oblate 
minima. For J = 2, 4 and 6 the lowest lying solutions {a = 1) are well lo- 
calized inside the prolate wells whereas the excited solutions (ct = 2) show 
collective wave functions reminiscent of a /3 vibrational state. On the right 
hand side of Fig. 7 wc have plotted the AMPGCM energies for cr = 1 and 
2 along with the / = projected energy. The AMPGCM energies have been 
placed along the q2o axis according to their "intrinsic average quadrupole mo- 
ment" which is defined as the average quadrupole moment weighted with the 
"collective wave functions" 

(920)a = y^rf920a20|5''"(920)r. (36) 

As a result of the mixing with the oblate configurations, the 0^ state has 
a lower intrinsic quadrupole deformation than the minimum of the J = 
projected energy but remains strongly deformed indicating that the N=20 
shell closure is broken in this nucleus. The intrinsic quadrupole deformation 
of the 2'1 , and 6^ states remains similar the the one of the corresponding 
minima of the AMP energy curves as a consequence of the localization inside 
the prolate wells of their their collective wave functions. The , ^J, etc, 
states have intrinsic quadrupole deformations similar to the i'l , 6f ones in 
good agreement with its quadrupole vibrational character. Another interest- 
ing result is the low excitation energy of the O2 state which is related to the 
strong prolate-oblate mixing. 

The results we have obtained with the AMPGCM for the 2+ excitation 
energy and B{E2,0f 2f) transition probability are 1.46 MeV and 395 
fm'', respectively. These results have to be compared with the angular mo- 
mentum projected results obtained with the approximate PBV calculation 
discussed in Sect. 3.5 (0.87 MeV and 593 fm^) and the experiment, 0.88 
MeV and 454± 78 fm'* . As a consequence of configuration mixing the 2"'" ex- 
citation energy increases substantially as compared with the AMP result and 
differs considerably from the experimental number. However, the AMPGCM 
B{E2) transition probability gets reduced with respect to the AMP result 
and gets closer to the experimental value. 

^^Mg though an spherical nucleus in the mean field approximation rep- 
resents the class of nucleus with a second minimum or a shoulder at an 
excitation energy comparable with the energy gain by AMP. A different case 
is provided by the nucleus ^'^^Pb already discussed or **Ca, see left panel of 
Fig. 8, which is also spherical at the mean field level but without any shoul- 
der at relevant. The effect of the AMP in this case is to produce two almost 
degenerated prolate and oblate minima at small q values. The behavior of 
(Jy) determine in part the energy gain of the ground state in the AMP, see 
Fig. 5 and the corresponding discussion, which is of the order of 3 MeV for 
^"^^Pb and of 1 MeV for ^^Ca. In this case a configuration mixing calcula- 
tion is called for which produces a ground state wave function which is a 
linear combination with similar weights of both minima and on the average 
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Fig. 8. The PES results for ''*Ca and ""Mg for 7 = 0, 2, 4, 6 and 8, along with the 
HFB energy (dashed curves) as well as the AMPGCM results (see comments on 
the right panel of Fig. 7 for an explanation) 



the ground state becomes again spherical, see Fig. 8. An additional energy 
gain (around 0.5 MeV in the case of ^^Ca) is also provided. For angular mo- 
mentum 2, 4, 6 and 8 we find in Fig. 8 that the prolate and oblate minima 
of the PES of ''^Ca take place at larger deformations with increasing spins. 
This produces larger barriers between the minima and as a consequence the 
mixing get smaller with growing spins and the wave functions remain well 
locaHzed in the corresponding minima. To complete our discussion we have 
plotted the mean field results, the PES and the AMPGCM energy levels in 
the right panel for the well deformed nucleus *"Mg, though the REC's are 
similar for both plotted nuclei, the effect of the projection on the mean field 
energy is quite different. The MF prolate ground state and the secondary 
oblate minima, after the AMP remain approximately at the same (J20 values 
and the same relative energy. The energy gain of the minima amounts to 3.5 
MeV. The structure of the MFA energy surface is also roughly maintained by 
the PES at higher spins. The configuration mixing calculations provide differ- 
ent results than in the spherical nucleus, one obtain two rather well defined 
rotational bands, i.e., with rather constant intrinsic quadrupole moment, and 
less mixing in the wave functions. The energy gain of the ground state by the 
configuration mixing amounts about 600 keV. 

In the following we will present results obtained with the Skyrme and the 
Gogny interaction. The calculations with the Skyrme interaction have been 
carried out mostly with the SLy4 [20] or the SLy6 [21] parameterizations of 
the force. In the AMPGCM ansatz a projection onto good particle number is 
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also performed. In the case of the Gogny interaction the parameterization DIS 
[22] has been used. A simultaneous particle number projection has not been 
carried out because the finite range of the force considerably increases the 
computational burden as compared to the Skyrme interaction. As the field 
is pretty new there arc not too many results and in some cases the nuclei 
computed with the Skyrme interaction are the same as the ones computed 
with the Gogny force. In those situations we will concentrate mainly in the 
Gogny results as the Skyrme ones are mostly in qualitative (and in some cases, 
quantitative) agreement with the former interaction. Most of the calculations 
have been focusing in light nuclei either with or without an excess of neutrons. 
Also calculations in heavier nuclei arc available in the neutron deficient Lead 
isotopes although some more results will come in the near future. 

We will start first commenting on the Skyrme results: the well deformed 
^''Mg nucleus was discussed in [48] with the SLy4 force and zero range pairing 
interactions (depending and not depending on the density) and a reasonable 
agreement was obtained for the ground state rotational band although the 
band was too stretched. The intra-band B{E2) transition probabilities were 
well reproduced. In [49] the doubly magic ^^O was studied with the SLy4 
parameterization and focusing on the deformed 0"*" excited states and the 
corresponding rotational bands built on top of them. A rather good agree- 
ment was found between the calculations and the experimental results. The 
structure of the first excited 0+ state was discussed in terms of multi particle- 
hole components and it was found that the structure of this state was a 4p-4h 
excitation in agreement with previous shell model explanations. It is inter- 
esting to point out that the ground state energy was lowered by 2.3 MeV 
with respect to the mean field as a consequence of the correlations induced 
by the AMPGCM method. The neutron deficient ^^^Pb, which is a very nice 
example of shape coexistence with and spherical ground state and two excited 
states below 1 MeV of prolate and oblate character respectively, was studied 
in [50] with the SLy6 parameterization and zero range density dependent 
pairing interaction. The excitation energy of the prolate 0+excited state was 
very well reproduced but this was not the case for the energy of the nearby 
oblate 0+state. The rotational band built on top of the prolate 0+ excited 
state came up too spread out as compared with the experimental data. The 
energies of both 0+ states are pushed up to higher energies (1.05 MeV and 
1.39 MeV for the prolate and oblate states, respectively) when the SLy4 force 
is used worsening the agreement with experiment. Finally, in [51] the nuclei 
^^S, ^^Ar, ^^Ar and "^"Ca were studied. The interest was here the normal and 
super-deformed bands recently found in those nuclei. The normal deformed 
and super-deformed band heads come up at a reasonable excitation energy 
in all the nuclei but the rotational bands built on top of them come always 
too spread out as compared with the experiment. However, the B{E2) values 
are in rather good agreement with experiment. The existence of SD bands in 
those nuclei is connected to a partial occupation of the /7/2 sub-shell. 
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The calculations with the Gogny force have focused in neutron rich Hght 
nuclei in order to investigate the erosion of the iV = 20 and TV = 28 shell clo- 
sures. To study the A'' = 20 region we computed [52] the Angular Momentum 
Projected energy curves for ^'^~^*Mg and 32-38gj g^j^^j from the minima of the 
projected curves wc determined the lowest state for each angular momentum. 
In this way we were able to demonstrate that the ground state of the magic 
^^Mg was deformed and traced back the breaking of the shell closure to the 
partial occupancy of the /7/2 intruder orbital. We also obtained a deformed 
ground state in '^'^Mg and for the other nuclei we concluded that they were 
showing shape coexistence. Concerning the excitation energy of the 2+states 
wc got a rather nice agreement with experiment as well as for the B(E2) 
transition probabilities. From our prehminary analysis it was clear that a 
proper description of those nuclei should involve configuration mixing in the 
context of the AMPGCM. This calculation was carried out for ^"-^Hlg [53] 
and for 32-36gj jg^j_ general rule, the inclusion of configuration mix- 

ing increases the energy of the 2+states as compared with the AMP results 
worsening the agreement with experiment but the B(E2) get usually reduced 
improving the agreement with experiment. It was estabhshed the spherical 
character of ^°Mg and the deformed one of ^^~^^Mg as weh as the spherical 
character of the ground states of all the Si isotopes studied. 

To analyze the erosion of the N = 28 shell closure we have considered [55] 
the isotones from ^"^Mg up to ^^Ca as weh as the Sulfur isotopes ^^~^^S. It is 
found that the drip line nucleus '*"Mg is prolate deformed breaking thereby 
the N = 28 shell closure, the nucleus ''^Si is also deformed but in this case in 
the oblate side. The nuclei **S and ^^Ar are found to show shape coexistence 
whereas ^^Ca is found to be spherical as expected. Concerning the Sulfur iso- 
topes they are found to be to a greater or lesser extent prolate deformed. The 
excitation energies of the 2+states and specially their transition probabilities 
to the ground state are rather well reproduced (the energies always come up 
rather high as compared with the experiment). The two neutron and two 
protons separation energies are also well reproduced and it is found that the 
effect of AMPGCM is only significant in a few cases while in the other the 
results are pretty close to the HFB ones. 

The Magnesium isotopic chain from A=20 until the neutron drip line 
A=40 was studied in [36]. The results obtained for the 2+, 4+ and O2 excita- 
tion energies as well as the B{E2,0'^ ^2+) are represented in Fig. 9 for the 
whole chain and compared with the experimental data. We observe that the 
isotopic trend is quite well reproduced but the 2+ and 4+ excitation energies 
come up too high whereas the Oj come up too low. Interestingly we are able 
to reproduce the shding down of the O2 excitation energy in going from ^^Mg 
to ^^Mg and due to the fact that whereas ^"'Mg is a well deformed nucleus 
^^Mg is a shape coexistent one [56]. One of the main results of our calcula- 
tion is that the nuclei from ^^Mg to *'^Mg show a prolate deformed ground 
state. To finish the discussion of the Figure let us finally comment that the 




Fig. 9. Comparison of the theoretical results (full line) for the 2+, 4+, and 
excitation energies and B{E2,0'^ — » 2"*") transition probabilities for the Mg isotopic 
chain with the experimental data (dashed line) 

experimental B{E2) transition probabilities are pretty well reproduced by 
our calculations. 

The AMPGCM method not only provides the excitation energy of the 
states but also the ground state energy and therefore one can compute two 
neutron or two proton separation energies and compare them with the HFB 
results. This comparison is made in Fig. 10 for the Mg [36] and Ne [57] iso- 
topic chains and there we observe how the AMPGCM method improves the 
agreement with experiment as compared with the HFB results. It is partic- 
ularly interesting the N = 22 result that is strongly related to the erosion 
of the A'' = 20 shell closure, the inclusion of correlations that made ^^Mg 
deformed are the responsible for the better agreement of the iS'2jv separation 
energy in this case. Finally, let us mention that according to our calculations 
the neutron drip line for the Mg isotopes is located at A^ = 28. For the Ne 
isotopes the AMPGCM predicts (contrary to the mean field) an stable ^^Ne 
in good agreement with experiment. 

Calculations similar to the ones presented for the Mg isotopes have also 
been carried out for the Ne isotopic chain in [57]. The conclusions are similar 
to the ones found for the Mg chain and they will not be discussed here. 

Finally, let us comment on the fact that independently of the force used 
the rotational bands obtained with the AMPGCM method are too stretched 
as compared with the experimental data. On the other hand the B{E2) tran- 
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Fig. 10. Two neutron separation energies S-zn as a funrtion of neutron number 
for the isotopic Mg and Ne chains. Results for the mean field and the AMPGCM 
calculations are present along with the experimental data 

sition probabilities compare pretty well with experiment. This "feature" of the 
method might have to do with the lack in the present formalism of additional 
correlations hke triaxialities, two-quasiparticle mixing and so on. 

6 Outlook and Future Developments 

In this paper we have described how to carry out Angular Momentum Pro- 
jection (AMP) with effective interactions of the Skyrme and Gogny type. We 
have illustrated the procedure with several examples. The main outcome of 
those calculations is that AMP strongly modifies the energy landscape as 
a function of the mass quadrupole moment of the system and therefore the 
naive procedure of Projection After Variation (PAV) (i.e. projecting after the 
intrinsic state is determined at the mean field level by optimizing the intrin- 
sic energy) leads in many cases to the wrong answer. On the other hand, 
the restricted Projection Before Variation (PBV) (i.e. the intrinsic state is 
determined by optimizing the projected energy) where only the quadrupole 
moment is allowed to vary yields much more consistent results. We have 
also learned that the projected energy landscapes have in many cases coex- 
isting minima and therefore it is important to consider configuration mixing 
which is implemented in the framework of the Angular Momentum Projected 
Generator Coordinate Method (AMPGCM). The results shown indicate the 
relevance of configuration mixing in many cases. We have also discussed ap- 
proximate ways to compute the projected quantities and in this context we 
have compared the present methodology with the results of the cranking 
model. All the methods considered have been illustrated with relevant ex- 
amples and at the end we have given a rather exhaustive account of the 
available theoretical results. As the field is rather new there are still further 
improvements to be made: a) It has to be explored if there are other relevant 
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degrees of freedom apart from the quadrupole moment that could be used 
in the context of the restricted PBV. b) The suitability of the method for 
the description of odd nuclei has to be explored, c) The octupole degree of 
freedom has to be incorporated in order to extend the present method into 
the realm of negative parity states, d) Triaxial and/or time reversal break- 
ing admixtures have to be incorporated, e) Approximations to compute the 
projected quantities have to be explored and their suitability assessed, and a 
long etcetera ! 
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